Electronic thermal conductivity as derived by density functional theory 
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Reliable evaluation of phonon thermal conductivity is of importance for understanding mecha- 
nisms of phonon scatterings and therefore further benefits optimizing figure-of-merit of thermoelec- 
tric materials. Usually, when experimentally deriving the phonon mediated thermal conductivity 
ftph = k — K e i from the measured total thermal conductivity k the constant Lorenz number Lo is 
chosen for estimating K e i according to the Wiedemann-Franz law for simple metals K e i = TLoo\ 
Present study demonstrates that such a procedure is not reliable when the Seebeck coefficient S 
becomes large which is just the case for a thermoelectric material. For a more reliable estimate it 
is proposed to use the modified Lorenz number Lo— S 2 which can be directly derived from the mea- 
sured Seebeck coefficient. Calculations by combining density functional theory with Boltzmann's 
semi-classical transport theory have been made for the clathrate type-I compound Ba8Au6-a;Ge4o+a; 
with x = and x = 0.27 corresponding to an electron doping of 0.8 states. For x = 0.27 the calcu- 
lated temperature dependent Seebeck coefficient agrees well with recent experiments corroborating 
the validity of the density functional theory approach. 

PACS numbers: 72.15.Jf, 72.15.Eb, 71.20.-b 
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Thermal conductivity plays an important role for the 
thermoelectric performance of a material as expressed by 
the figure-of-merit ZT = TS 2 a / (n e i+K p h) which includes 
the Seebeck coefficient S, the electrical conductivity a, 
and the thermal conductivity k = n e i + n p h summing 
up the contributions of electronic states and phonon me- 
diated processes. Consequently, a low thermal conduc- 
tivity in combination with large values of S and a are 
desirable in order to achieve large values of ZT. Consid- 
erable efforts for lowering k by reducing n p h were made by 
utilizing structural properties, such as building up super- 
lattices and incorporating suitable filler atom into 
structural cages [5rJl2j . These concepts rely on the strong 
scattering of heat-transporting phonon modes. Often the 
effect of phonon scattering is taken into account by es- 
timating K p h with phenomenological formulas. However, 
neither K e i nor n p h are directly measured. Rather, k p u 
is extracted by subtracting n e i from the measured total 
thermal conductivity, i.e., n p h ~ K exp — n e i in which the 
electronic thermal conductivity is derived by employing 
the Wiedemann-Franz law for simple metals, n e i « TLqct 
[fj|-[lll Il3rl2l| . In this expression, Lo is a universal con- 
stant and does not depend on temperature and on the 
properties of the material under investigation. As will be 
shown by a density functional theory study for a ther- 
moelectric materials this assumption of a simple metal 
may lead to unreliable estimates in particular when the 
Seebeck coefficient of the material is large. 

The present theoretical study of the electronic ther- 
mal conductivity is based on the same density functional 
theory (DFT) concept as used for first-principles calcula- 
tions of Seebeck coefficients (for example, see Ref. [22} ) ■ 
The Wiedemann-Franz law is generalized by introducing 
a material and temperature dependent Lorenz tensor L. 



For its calculation Boltzmann's transport theory is ap- 
plied in combination with electronic properties as derived 
by a DFT calculation. This procedure is applied for the 
clathrate type-I compound BagAue-aGe^+a, which is a 
prototypical thermoelectric material and for which also 
recent measurements of Seebeck coefficients are available. 
The Wiedemann-Franz law is generalized to 



K e i = TLcr 



(1) 



in which the tensor n e i is related to the conductivity ten- 
sor a of Eq. IA.2I via the Lorenz tensor L. These quanti- 
ties as well as the Seebeck tensor S (Eq. I A.3|) are second 
rank tensors [23j . In accordance with Boltzmann's trans- 
port theory one derives [24| 



«ei = i (K 2 - K^KcT 1 ) 



(2) 



for which Eq. IA.2l is utilized for the definition of the ten- 
sor K n . The Lorenz tensor L can now be formulated 

as 



Li 



whereby Li is expressed as [24 



Li 



1 



e 2 T 



rK 2 K 



(3) 



(4) 



For a free-electron like metal the the second term at 
the right hand side of Eq. [3] is negligible. In this case the 
Lorenz tensor can be approximated by L w Li. Inspect- 
ing Eq. [3] it becomes obvious that the deviation from 
the free-electron like behavior is given by — S 2 . At low 
temperatures for a simple metal the original Wiedemann- 
Franz law is a reasonable approximation, i.e. Li ~ Lq, 
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whereby Lo would be a tensor with constant coefficients 
L = = 2.44 x lQ- 8 Wn/K 2 . 

When assuming the constant relaxation time approxi- 
mation -as it is also standard for first-principles calcula- 
tions of the Seebeck coefficients [22|]- t cancels out in the 
components of L since it appears in the numerator and 
denominator according to Eqs|3] and SJ This is also the 
case for the Seebeck coefficients given by Eq. IA.31 For 
the following discussion it should be noted that crystal 
structure of the material under study is of cubic sym- 
metry. As a consequence of this high symmetry all the 
second rank tensors as mentioned above are diagonal and 
the three diagonal coefficients are equal. Therefore, only 
one coefficients needs to be considered for each tensor. 
However, the derivations and calculations can be done 
for a general crystal symmetry and for tensors with less 
symmetry and more components. If the tensors are sym- 
metry averaged again only one coefficients needs to be 
considered. 

What remains for the calculation of the integrals ap- 
pearing in Eq. IA.2I is the calculation of the electronic 
structure and related properties of the actual solid ma- 
terial, which was done for Ba 8 AueGe4o by DFT. For 
that, the Vienna Ab initio Simulation Package (VASP) 
[25l |26| was used for which the pseudopotentials were 
constructed according to the projector augmented wave 
method [27!, 28] . The exchange-correlation functional was 
parametrized in terms of the local density approximation 
according to Ceperley and Alder [29]. The valence state 
configuration for the construction of the pseudopoten- 
tials included the 5s, 5p and 6s states for Ba, the 6s and 
5d states for Au, and the 3d, 4s and 4p states for Ge. 
For the Brillouin zone integration a 5 x 5 x 5 grid of k- 
points was chosen, which results in well-converged total 
energies and fully relaxed structural parameters as given 
in the caption of Fig. HJa). Electronic transport proper- 
ties were derived by solving the semi-classical Boltzmann 
transport equation within the constant relaxation time 
approximation as imp lemented in the program package 
BoltzTrap [H, [3(| HJ . For that purpose the Kohn-Sham 
eigenvalues £j(k) were determined on a very dense 25 x 
25 x 25 grid of k-points. 

Fig. Q^a) shows the density of states (DOS) around 
Fermi energy. A gap of about 0.3 eV occurs about 0.15 
eV above the Fermi energy Ep for the undoped case. The 
electron doping of AA = 0.8 places E^ very close to the 
gap where the DOS changes rather strongly. For such 
a situation one would expect a large Seebeck coeffcient, 
as is indeed the case (see inset in Fig. [5]). By adding 
AN = 2 electrons the corresponding Ep falls into the 
gap. For T=300K FigQJb) depicts the coefficients L and 
Li of the cubic Lorenz tensors according to Eqs. [3] and 
[3J as functions of the doping. It is obvious that the de- 
viation between the approximation Li and L of the full 
calculation increases strongly the closer the doping is to 
AN = 2. It should be noted that for AN = 2 the co- 
efficients L and Lj are undefined because no electronic 
states are available for transport and, subsequently, the 
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FIG. 1. DFT results for BagAu6Ge4o: (a) density of states 
(DOS) vs. energy E whereby the Fermi energy for the un- 
doped case corresponds to E = 0. The Fermi levels for the 
undoped (AN = 0.0) and an electron doped (AN = 0.8 or 
x — 0.8/3) compound are indicated by a full and dashed 
line, respectively. Doping was modelled by the rigid band 
approximation, i.e. by a suitable shift of the Fermi energy 
with respect to the undoped case but leaving the DOS un- 
changed otherwise. BasAueGe-io has the type-I clathrate 
structure(cubic unit cell, space group Pm-3n). The relaxed 
structural parameters are 10.70 A for the cubic lattice pa- 
rameter and the site specific coordinates are 2a (0,0,0) and 
6c (0.25,0.5,0) for Ba, 6d (0.25,0.5,0) for Cu, 16i (0.183, 0.183, 
0.183) and 24k (0, 0.117, 0.309) for Ge. (b) Coefficient L of 
the cubic Lorenz tensor and its approximation Li (Eqs. [3] HJ) 
at 300K vs. rigid-band doping AN. Negative(positive) values 
correspond to hole(electron) doping. For AiV = 2 the Fermi 
energy falls into the gap of the DOS. The constant Lorenz 
number Lo of the Wiedemann-Franz law is indicated by a 
dashed horizontal line. 



electrical conductivity a is zero at lower temperatures. 
The comparison between Li and L demonstrates the ef- 
fect of the contribution of the Seebeck coefficient, -S 2 . 
Significant deviations occur in the range of < AA~ < 4. 
For dopings outside this critical region the coefficients 
L and Li coincide and their values are comparable to 
the constant Lorenz number but still Lq might be off by 
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FIG. 2. DFT results for BagAue-xGe-io+z: Coefficient of 
the Lorenz tensor as a function of temperature for no doping 
(black solid line) and for a doping of AiV = 0.8 electrons (blue 
dashed line), see Fig. [T] The constant Lorenz number Lo of 
the Wiedemann-Franz law is indicated by a dashed horizontal 
line. The inset shows the corresponding DFT-derived Seebeck 
coefficients in comparison to the experimental values of Ref. 
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FIG. 3. DFT results for BagAus-xGe^+i with a doping of 
AiV = 0.8 electrons: electronic thermal conductivities K e i as 
calculated by using the coefficient L of the full Lorenz tensor, 
the coefficient Li of its free-electron approximation, the con- 
stant coefficient Lo and the suggested approximation Lq-S 2 . 



20%. When the Fermi level is sufficiently far away from 
the gap, L has almost the same value as Li. The results 
suggest that the approximation of Eq. (0} is not valid for 
thermoelectric materials with large Seebeck coefficients 
and it will lead to an overestimation which might be 
large. 

Fig. [3] shows the temperature dependent coefficients 
Li and L for no doping and AiV = 0.8. The inset 
compares DFT derived and experimental Seebeck coef- 
ficients revealing rather good agreement between calcu- 
lation and measurement for the DFT data corresponding 



to AiV = 0.8. This change in valence electron num- 
ber refers to the slightly off-stoichiometric composition 
BasAu6-a;Ge4o+a; with x = 0.8/3, which is within the 
error bar of the measured stoichiometry [32j. Both Li 
and L approach the Wiedemann-Franz limit Lo at low 
temperatures but deviate significantly at elevated tem- 
peratures. Overall, L and Li exhibit a strong tempera- 
ture dependency and the deviation of the constant value 
L gets large in particular around 700K, which is in the 
range of temperatures as they occurin technological ap- 
plications. 

For finally deriving the electronic thermal conductiv- 
ity K e i from Eq. Q] the full Lorenz coefficient L together 
with its approximations Li and Lo is used. The results 
of in Fig. [3] reflect the behavior of the Lorenz coeffi- 
cients in Fig. [2] showing again very significant differences 
up to 40% between the full calculation and the result 
involving the approximations Li and Lo. Clearly, just 
using the classical Wiedemann-Franz law (i.e. Lo) may 
lead to rather unreliable values for the electronic thermal 
conductivity. If then K e i is used for deriving the lattice 
thermal conductivity from k = n e i + K p h the results for 
K p h will also be unreliable. Comparing the rather simi- 
lar behavior of K e i in Fig. [3] as calculated with Li and 
L it seems a suggestive strategy to correct the Lorenz 
constant Lo by subtracting S 2 as is done in case of Li in 
Eq. [3] Doing this gives good agreement with results as 
derived from the coefficient L of the full Lorenz tensor, 
as shown in Fig. [3] Therefore, if the temperature depen- 
dent Seebeck coefficient S(T) of a material is measured 
in combination with the electric conductivity by applying 
the corrected Lorenz number Lo— S 2 a more reasonable 
estimate is achieved for the electronic thermal conduc- 
tivity and, consequently, for the lattice thermal conduc- 
tivity than it is the case when only Lo is used. Achieving 
reasonable estimates for n p h by reasonable estimates of 
K e i is °f importance for optimizing transport properties 
of thermoelectric materials. 

In summary, in this work separating the electrical and 
lattice thermal conductivities, i.e., n e i and K p h, using 
Lorenz number is studied. Calculations were performed 
by hybridizing DFT and Boltzmann transport theory, 
which shows that the Lorenz number in general is a ma- 
terial and temperature dependent quantity rather then a 
constant value as it is usually derived from the metallic 
limit. In particular, significant deviations from the con- 
stant value are yielded when compound is a good ther- 
moelectric material. Using the constant Lorenz number 
causes severe discrepancies in the resultant K e i from that 
using fully calculated Lorenz number, which become sig- 
nificant at higher temperatures. In the present study, an 
error about 40% at 1000 K was obtained for BasAu6Ge4o 
. The results suggest that one has to rethink the way of 
extracting K e i and K p h from experiment in the usual way. 
To remedy this situation, we propose that instead simply 
using the constant Lorenz number for all materials over 
all temperatures, one needs to subtract the square of See- 
beck coefficient, an element quantity from thermoelectric 
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measurements, from the constant, which can give rise to 
a more precise evaluation of the lattice thermal conduc- 
tivity K ph . 
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Appendix 

The relation needed in Boltzmann's transport theory 

is 



K » = 4^3 12 r *( k ) v *( k ) ® v i( k )M k ) ~ /")" (A.I) 

in which Tj(k) is the relaxation time of the electronic 
states with band index i, energy eigenvalue £j and band 
velocities for wave vector k, whereas /(/i, T denotes 
the Fermi-Dirac distribution function for the chemical 
potential [i at temperature T. 

The electrical conductivity a is given by 

a = e 2 K (A.2) 
and the Seebeck tensor S is defined by the relation. 

S = 4;KiK - 1 . (A.3) 
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